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Abstract 

The quantum Fourier transform (QFT) plays an important role in many known quantum algorithms such 
as Shor's algorithm for prime factorisation. In this paper we show that the QFT algorithm can, on a 
restricted set of input states, be de-quantised into a classical algorithm which is both more efficient and 
simpler than the quantum algorithm. By working directly with the algorithm instead of the circuit, we 
develop a simple classical version of the quantum basis-state algorithm. We formulate conditions for a 
separable state to remain separable after the QFT is performed, and use these conditions to extend the 
de-quantised algorithm to work on all such states without loss of efficiency. Our technique highlights the 
linearity of quantum mechanics as the fundamental feature accounting for the difference between quantum 
and de-quantised algorithms, and that it is this linearity which makes the QFT such a useful tool in quantum 
computation. 
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1. Introduction 

The quantum Fourier transform ( QFT) plays an important role in a large number of known algorithms 
for quantum computers [1]. It plays a central role in Shor's algorithm for prime factorisation [2| and is 
often thought to be at the heart of many quantum algorithms which are faster than any known classical 
counterpart. However, following on from recent results relating to classical features of the QFT algorithm 0- 

we will argue that the QFT algorithm itself is classical in nature. 

The process of de-quantising quantum algorithms into equivalent classical algorithms is a powerful tool 
for investigating the nature of quantum algorithms and computation. Few general results are known about 
when such de-quantisations are possible and the power of quantum computation compared to classical 
computation. In this paper we show how the QFT algorithm can be de-quantised into a simpler, more 
efficient, classical algorithm when operating on a range of input states. While the de-quantised algorithms 
themselves are of interest, they also allow us to gain insight into the nature of the QFT. We will argue that 
it is the linearity inherit in the unitary quantum computational model which makes the QFT such a useful 
tool, rather than the nature of the QFT itself. 

In Section[2]of this paper we overview the basic QFT theory and present the QFT algorithm in a compact 
form which allows us to move away from the restrictions imposed by the circuit layout. In Section [3] we 
overview the de-quantisation procedure and de-quantise the QFT algorithm acting on a basis-state input. 
In Section 2] we explore the entangling power of the QFT and determine conditions for when a separable 
input state remains unentangled by the QFT, before presenting a de-quantised algorithm that works on such 
product-state inputs. In Section [5] we discuss why de-quantisation of the QFT is possible and note some 
common misunderstandings about the QFT which contribute to this. 
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2. Discrete and Quantum Fourier Transforms 



The discrete Fourier transform (DFT) on which the QFT is based is a transformation on a g-dimensional 
complex vector x — (/(O), /(I), • ■ • , /('?~ 1)) into its Fourier representation x — (/(O), /(I), • ■ • , f{l~^)) Q- 

/(c) = ^|]e2— (1) 

for c G {0, 1, . . . ,q — 1}. The QFT is similarly defined so that the transformation acts on a state vector in 
g-dimensional Hilbert space, Hq. In quantum computation we work with a state vector defining a register 
comprising of n two-state qubits, so we will only consider the case that q = 2" from this point onwards. We 
will use the convention that n is the number of qubits while iV = 2" is the dimension of Hilbert space the 
n qubits are in. This means that the QFT, denoted Fq, acts on the N amplitudes of a particular n-qubit 
state, i.e. 

N-l N-l 

5]/(a)|a)^5]/(c)|c). (2) 

a=0 c=0 

The QFT hence transforms a state so as to perform a DFT on its state vector. 

As a result of the linearity of quantum mechanics, in order to compute the QFT we only need to 
design an algorithm to transform a single component of the state vector. This is because an arbitrary state 
IV'Tv) = '^a=o /('^) 1*^) transforms as: 

N-l N~1N-1 N-l 

Fn IViv) = E /(«)^^ l«) - 7^ E E Ic) = E /(c) 1^) " 

a=0 ^ a=0 c=0 c=0 

Hence we arrive at the standard definition of the QFT as the mapping Q 

N-l 



a 



Fn. 1 



2iTiac/N 



N ^ 

c=0 



e'"'"'!''' |c) , (3) 



with a G {0, 1, . . . , — 1}. Keeping in mind that we are dealing with registers composing of qubits, we can 
decompose a (and similarly c) into its binary representation so that a — 2"^^ai + 2"^^a2 + - • • + 2^a„_i + 2°a„ 
and \a) — \a1a2 ■ ■ ■ an)- By denoting a = 0102 • • • On and a/2" — Q.aia2 • • • we observe that 

g27riac/2" _ g27ria(2""^ci+2""^C2H |-2"c„)/2" 

^27rz(aia2'"aii)ci/2^g27rz(aia2'"aTi)c2/2^ ^TTi{a\a2 ■■■an)cn/2^ 

_ g27ri(ai---a„_i.a„)cig27rj(ai---a„_2-a7i-ian)c2 . . . g27ri(0.aia2---a„)c„ 

Noting that for any decimal x.y we have e^'^'>-i^-v) = (•g27ri^Xg2Tri(o.a) _ ^2-Kt{o.y) ^ ^]^,^^ ^^y^ |.j^g fractional 

part of (ai • • • a„_j.a„_j+i • • • an)cj is of any significance in the exponent of QII] Hence, we find 

Using this decomposition we can write as a product state of individual qubits, 

N-l 

e^'^Wa" |c) ^ (|o) + e^'^^f"-"") |1)) • • • (|0) + e^^^^O-'^i-'^") |1)). (5) 

c=0 



^This technique of removing factors of (e^'^')'' for fe £ N will be commonly used throughout this paper to reduce formulae. 
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1^2) 



H - R2 



Rn-l - Rn 



\an^i 
la. 



n/ 



H 



Rn-1 



Rn 



- H - R: 



|0) + e27ri(0.aia2---a„) j^^ 
|0) + g27ri(0.a2a3---a„) j^^ 

|0) + e27ri(0.a„_ia„) jj^^ 



H - |0> +e2"(0-''") |1) 



Figure 1: The standard quantum circuit for the QFT. The output normalisation factors of l/\/2 and swap gates to reverse 
qubit order are omitted. 



The quantum algorithm to implement the QFT follows directly from this factorisation. The circuit for the 
algorithm is shown in Figure [TJ The algorithm can be written explicitly as follows 0] : 



Quantum Fourier Transform 

Input: The state \a) = \ai) \a2) ■ ■ ■ |a„). 

Output: The transformed state -;^(|0) + e^'^'^O-"") |1)) 

1. For J = 1 to ri, transform qubit \aj) as follows: 

H 



(|0) 



g27ri(0.ai 



2. 
3. 
4. 



For fc = j + 1 to n: 

^(|0) + e2'^*(0-°^-°fe-i) |1)) ^ -i=(|0) + e2'^»(0-a.-a/c-i"fc) |i)) where Rk 
is the unitary /c-controlled phase shift: 



Rk — 



/I 

10 

1 

VO e2-V2' 



End For. 

Reverse the order of the qubits. 



7. End For. 

Clearly this produces the state ([5]) and requires O(n^) unitary Rk or H gates to run. 

There are a few important notes about the QFT which should be made. While both the DFT and 
the QFT act on vectors in a complex vector space, the DFT acts on an abstract, mathematical vector, 
whereas the QFT acts on a physical state which we mathematically represent by a vector in "Hat. The 
subtle difference here is that with the classical DFT, we can read the values of all 2" Fourier coefficients 
/(c) by simple inspection of the transformed vector. With the QFT, the resulting state embeds all 2" 
coefficients as amplitudes for the 2" states of an n-qubit system. However, the collapse of the superposition 
upon measurement means that it is impossible to measure the amplitudes of a quantum state without an 
ensemble of such states to make a statistical approximation of the amplitudes from 0] , and detecting phase 
differences between states is even more difficult. Hence, the quantum state ^ contains all the information 
of the classically transformed vector, but it is inaccessible to measurement. The main use of the QFT is then 
as a tool to extract information embedded in the relative amplitudes of states as opposed to determining 
the coefficients themselves. 

Another result of this is that the efficiency of the QFT (O(n^) as opposed to the DFT which is 0(n2")) 
is in some sense due to the ability to perform the transformation and utilise the information in the phases 
without measuring the state. Evidently, any algorithm requiring measurement needs exponential time 
(there are 2" coefficients to measure), so even if quantum mechanics would allow us to measure the Fourier 
coefficients in state ([2]), doing so would take 0(?i2") time: 2" coefficients, n qubits each. Making use of 
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this embedded information while avoiding measurement is certainly an important part of the fine art of 
developing algorithms in quantum computing. 

3. Initial De-quantisation Investigation 

Having presented the QFT, there are some issues to be brought to light. The decomposition of the 
transformed state ([3]) (shown in ([5])) is evidently not entangled, and the separability of the state would lead 
us to believe that the QFT algorithm producing it could be simulated efficiently in a classical manner , 
and there are certainly results towards this. 

It was realised shortly after the discovery of Shor's algorithm that the QFT could be computed in a 
semiclassical manner [5]. By using classical signals resulting from quantum measurements, one can perform 
the QFT on a state using classical logic and one-qubit gates (instead of the usual two-qubit controlled- 
phase-shifts). This method gives the same resulting probability distribution as the quantum algorithm, but 
destroys the state's superposition as it relies on irreversible measurements. As a result, this is only useful in 
an algorithm in which the QFT directly precedes measurement. Shor's algorithm happens to be of exactly 
this nature, but this is only an initial step towards true classical simulation. 

Much more recently, classical simulations of the QFT have been studied from the viewpoint of simulating 
the circuit in Figure [T] by exploiting the bubble- width of the quantum circuit (sl and the tensor contraction 
model 0. The bubble- width approach uses a slightly modified version of the QFT circuit which is of 
logarithmic bubble- width and simulates this circuit. The tensor-contraction model also focuses on the 
circuit topology, but relies on associating a tensor with each vertex in the circuit, then cleverly contracting 
the tensors into a single rank-one tensor. Both these methods work on separable input states, but are 
sampling based forms of de-quantisation flT| in the sense that a final measurement is assumed and an 
output is classically sampled from the correct (calculated) probability distribution. This makes these de- 
quantisations less general than might be desired and difficult to apply when the QFT is used, as it often is, 
as a part of a larger quantum algorithm. This is because in these cases measurement cannot be assumed 
after the QFT, and the de-quantisation must be cleverly and non-trivially composed with a de-quantisation 
of the rest of the algorithm to be applied. 

Working with the circuit topology, while beneficial for some purposes, also seems to overcomplicate 
matters and restrict generalisation when it comes to classical simulation. We will explore simulations of the 
QFT in a different light, more along the lines of the de-quantisation explored previously by Abbott Q and 
Calude fl3\ which aim to provide stronger (not sampling based) de-quantisations when possible. 

3.1. De-quantisation Overview 

The idea behind this de-quantisation procedure is that qubits which are separable exhibit only superpo- 
sition and interference. These properties are the result not of non-classical features of the qubits, but rather 
of the two-dimensionality of the qubits. By using classical, deterministic two-dimensional bits instead of 
qubits, the same behaviour can be exhibited without the difficulties imposed by measurement and prob- 
abilities. Not all algorithms fit within this paradigm, but there are many which can be tackled with this 
approach. Algorithms which use measurement as a fundamental part of their procedure are examples of those 
which are not so well suited, and sampling-based techniques are more suitable in these situations. Finding 
when these stronger de-quantisations are possible also gives insight into the power of particular quantum 
algorithms llj, as this reflects to some degree the amount the algorithm utilises the possible advantages of 
quantum mechanics. In cases where entanglement is bounded [l^ . we can use this de-quantisation procedure 
to produce classical algorithms which are as efficient as their quantum counterparts. This procedure was 
explicitly examined further [^, [l^] when applied to the Deutsch- Jozsa problem 13, 14\ , where complex num- 
bers were used as classical two-dimensional bits. In this paper we will apply this de-quantisation procedure 
to the QFT, but because the amplitudes we need to represent in the QFT algorithm are complex-valued, we 
cannot use complex numbers as our two-dimensional bits. There is no problem though with simply using 
two- valued vectors as our classical bits, so we will employ this procedure. 
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3.2. Basis-state De-quantisation 

The de-quantisation for a basis-state needs only to simulate the transformation defined in ([21). As a 
result of the decomposition in ([S]), the effect of the QFT on the jth qubit is easily seen to be 

|a,)^-L(|0)+e^-(°— (6) 

The difficulty in implementing this in a quantum computer is that the phase of a qubit needs to be altered 
depending on the values of the other qubits without altering them - that is why it is not helpful to express the 
quantum algorithm as we have done in ([S]) - and the circuit of controlled-phase-shifts is required to implement 
this. The information is spread over the input qubits and must be obtained without measurement. In the 
classical case there are no such restrictions on measurement, so de-quantisation should only require directly 
implementing (jS]). However, evaluating the complex phase for each of the n qubits takes 0{n) time, leading 
to a 0{n^) procedure. This can be reduced to 0(n) by calculating each phase dependent on the previous 
one. To do so, let uij be the jth phase factor and note the following: 

— g27ri(0.a„_j + i---a„) 

_ g27ri(0.a„_j + i)g27ri(0.a„_j+2---an)/2 



and 



uji = e2"("-''") = (-1)°", 



where by the square-root we mean the principal root. The square-root of a complex number such as ujj can 
be calculated independently of n. Specifically, if we have s + ti = \/b + di with the further requirement that 
for a root of unity + d'^ — 1, then [l^ : 

1 / r sgnfd) , 

V 2 V 2 

where sgn{d) — d/\d\ is the sign of d. The efficient de-quantised algorithm is then the following: 

Basis-state De-quantised QFT 
Input: The binary string a = 0102 . . . a„. 

Output: The n transformed two-component complex vectors ^162 . . . 6„. 

1. Let a; = 1 

2. For j = 1 to n: 

3. Set w = 



4. Set bj^^x 

5. End For 



This algorithm produces vectors mathematically identical to the state-vectors in ([3|) and ([5|) produced by 
the QFT, but is computed classically in 0(n) time - more efficient than the quantum solution and simpler 
too. This is primarily because the quantum circuit is constructed subject to the requirement of computing 
the QFT without any intermediate measurements. As a result, the quantum algorithm corresponding to the 
circuit must conform to this too, making it more complex than an equivalent classical algorithm need be. 

A classical algorithm has the further advantage over the quantum algorithm acting on a basis-state 
that measurement of the resulting state can be performed at will, and any required information is easily 
accessible. In the quantum algorithm only a single state can be measured, and no information about the 
amplitudes (and thus the Fourier coefficients) can be determined from a single QFT application. While 
this classical algorithm is no faster than the well known fast Fourier transform (FFT) for calculating all the 
coefficients, it may be advantageous if only some coefficients are required. 
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The ability to de-quantise the QFT acting on a basis state is not particularly surprising. This is equivalent 
to the classical DFT acting on a vector with only one non-zero component, producing a fairly trivial and 
easily computed output. However, this highlights a little more deeply some common misconceptions about 
the QFT. Because of the linear, unitary evolution of quantum mechanics, the action of the QFT on a basis 
state shown in ([3]) is often taken as the definition of the QFT. While this suffices as the definition for the 
purposes of the quantum algorithm, it is important not to forget that the actual definition of the QFT is 
that given in 1^. When considering classical simulations of the QFT this is even more important, as the 
action of the QFT on a basis state and the corresponding circuit no longer immediately allow us to compute 
the complete QFT; indeed it would take 2" iterations of a classical algorithm simulating the basis state 
behaviour to compute the complete QFT. 



4. Product-state De-quantisation 

Here we consider the possibility of extending the de-quantisation to work on a wider range of input 
states, resulting in a less trivial de-quantisation. If the input state is entangled then it is clear that the de- 
quantisation is not easily extended, as the method used for the basis-state algorithm relied on the separability 
of the input. In such a situation, any de-quantisation attempt would need to involve a different method and 
work directly from the QFT definition, 

It is not immediately clear that the basis-state de-quantisation, which is based on could not be 
extended to work on arbitrary separable input states. This idea is strengthened by the fact that we used 
the single-qubit formula ^ to perform the basis-state de-quantisation. However, this implicitly relies on 
the other qubits in the input state having a definite value, but in the general separable input case this is 
not necessarily the case. Indeed, the QFT is readily seen to entangle separable input states, e.g: 

10) = i= |o) (|o) + |i)) ^ (^100) + ^ 101) + ^ |ii; 

A de-quantisation for arbitrary separable input states is thus not possible in the same way as it was for basis 
states. However, we will investigate the entangling power of the QFT in order to determine the set of states 
which are not entangled by the QFT, and present a de-quantised algorithm which works for such states. 

4-.1. General Separability Conditions 

As in the entanglement investigation of the Deutsch- Jozsa problem Q , we will make use of the separability 
conditions for a qubit state presented in |16], although unlike the Deutsch- Jozsa problem our situation 
permits the possibility of states with zero-valued amplitudes, complicating the conditions somewhat. The 
key definitions and theorems we require to determine the separability of a state will be bricfiy presented, 
while [l^ should be consulted for proofs and discussion. 

Definition 1. The amplitude abstraction function A : Hn ^ {0, 1}^ is a function which, when applied to 
a state I^Aat) = J2i=o^ yields a bit string x — xqXi . . .x^^i such that for < * < A^ ^ 1, a^i = if 

Ci — and Xi = \ otherwise. 

Definition 2. The set Bn C {0, 1}^ of well-formed bit strings of length = 2" is defined recursively as 

62 = {01, 10, 11}, B2N = {O^x, xO^, xx\x& Bn}- 
Definition 3. The set of well-formed states is the set 

Vn - {\-^n) G Un I A{\iIjn)) e Bn}. 

Intuitively, a state is well-formed if the zero- valued amplitudes are distributed such that it is a candidate 
to be separable; if a state is not well-formed it is guaranteed to be entangled. In order to determine if a 
well-formed state is separable, we require two further definitions. 
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Definition 4. For each set of well- formed states Vn, there exists a family of zero deletion functions {Vk '■ 
Vjv — > Hk \ K = 2'',! < k < n}, such that for a well-formed state \iPn) — J2iLo^ \i) e Vat, I?A'(|V'Af)) — 
IV'k) — J2f=o^ ^'j b)j (Wk)) — Cj is the jth non-zero amplitude of IV'a')- 

Definition 5. A state \iPn) = ^iLa^ Ci \i) is pair product invariant if and only if for all j £ {2, . . . ,n} and 
all i G {0, . . . , J/2 — 1} CiCj^i-i — dj, where each dj is a constant and J — 2^ . 

As a concrete example to help understand pair product invariance, consider the cases of n = 2 and n = 3. 
For n = 2, \tp4) = X]?=o '^^ P^^'^ product invariant if the well known condition cqCs = ciC2 holds. For 
n = 3, \4>s) — X]I=o '^i I*)' require this same condition, cqCs = ciC2, as well as the further condition that 

CoCy = CiCg = C2C5 = C3C4, to hold. 

The following theorem from [16] can be used to determine if an arbitrary rt-qubit state is separable or 
not by checking the non-zero amplitudes of the state vector are pair product invariant. 

Tiieorem 1. Let \^n) be an n-qubit state for which the bit string AdipN)) contains K ones. Then \4'n) is 
separable if and only if H^n) G Vat and Vk HiPn)) is pair product invariant. 

4.2. QFT Separability Conditions 

We wish to consider the case that a separable n-qubit input state remains separable after the QFT has 
been applied to it. In order to do so, first let us consider the action of the QFT on the separable input state 

1^^) = ® (^2^ ® • • • ® (^2^ = (/(o), /(I), ...j{N-i)f. 

Note that each /(c) can be written as a product of amplitudes as /(c) = 0102 . . . a„, where each e {ai, /3i}. 
We will use the notation fj{c) to mean OjUj^i . . . a„, and thus /(c) — /i(c) ~ 01/2(0) etc. Because of the 
structure of the tensor product, for < j < n and c < 2""-', fj{c) — ajfj+i{c) and /, (2"^-' -}-c) = I3jfj^i{c). 
The amplitudes of the transformed state \^n) = (/(0),/(l), . ■ . ,f{N— 1))^ are given by (P), which can, for 
a separable input, be rewritten in the more useful form 

^ a— 

N/2-1 N/2-1 
= ^ai e'"^"'/'^/2(a)+/3l ^ ^2MN/2+a)c/N ^^^^^ 

N/2-1 

= (c,i+e--/3i) e2--/^/2(a) 

a=0 

= ^(ai + e"=/3i)(«2 + e-^^/2/32) • • • (a„ e"^/^"-^^^) 
V iV 

Tl 

^^\{{a,+e-'^-'p,). (7) 

This factorised form of the transformed Fourier coefficients allows us to determine conditions for when the 
transformed state is well-formed by giving restrictions on the distribution of zeros amongst the amplitudes, 
a result of the fact that /(c) = if and only if one of the factors in ([7]) is zero, and this is a significant step 
towards determining if a state is separable, and thus de-quantisable. 

Lemma 1. Let \iI}n) be a separable input state and \i^N) — -Fjv iV'iv) be the transformed state. Then the 
following three conditions are equivalent: 

(i) \iPn) G Vtv, *.e. the transformed state is well- formed. 
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(ii) There exists a k < n such that the set 

Cj ^{c\yi<j (a/ + e""/2'-'^^ ^0^1 

is non-empty for all 1 < j < k and empty for k < j < n. 

(iii) (30 < fc < n)(3ai . . .Ofe e {0, l}'^') (vi <j<k aj = ^p=i ''"/^'^''/Jj 

A(Vafc+i . . .a„ e {0,l}"-'^-)(Vn > j > k) \aj ^ e'^* ^p=i '^''/^^"''^^.j ^ 

Proof. Q For any x G Bn, Definition [5] ensures that the number of ones in x, = 2™ for 

some m < n, and hence the number of zeros, #o(a^) — 2" — 2™ = X^fcTi™ 2"^'. If |Cj| ^ then there exists 
a c' G such that c' < 2^ and /(c') = 0. But we must also have f{m2^ + c') = for < m < 2"--'' - 1 and 
hence \Cj \ — 2"^-'. Also note that each Cj is disjoint by construction, and /(c) = c G Cj for some j. 
For a well-formed state we thus require that for some m, 

n—m n 

#0 (^(I^at))) = 5] 2"-' = 5] |C,H E 2"^'' 
'=1 j=i iHCjMo 



which is satisfied if and only if Ci . . . are non-empty and Ck+i ■ ■ - Cn are empty, with k = n ~ m. 

(El) =^ 0: In the first K = 2'' amplitudes, 2''-''J2j<k = Ej=i 2''"^' = iiT - 1 of them are 
zero. Let /(c') be the single one of these non-zero amplitudes. Then, by symmetry, f{dK -t- c') 7^ for 
< d < 2"-'= - 1. Thus, AQiPn)) = x'^"'\ where x e {0,1}^ and #i(x) = 1. Any such x is clearly 
well-formed, and thus the state l-^Ar) is also well-formed. 

dn]) ■i=> ([m]): Note that J2p=i flp/2-'^^ = X]p=i o-p'^^^^, and we will proceed by induction for j < k. 
Since ai = e""i/3i ai e^^^+^i)/?! = 0, such an ai S {0,1} exists if and only if \Ci\ ^ 0. Now, 

assume that for aU 1 < / < j < A:, a; = ^p=^ ""^^ and |C/| 7^ 0. Then 

= ^p=i ""2"" ^ + e#^^^'"'+^^=i '^''^""'^/S^- = 0, 

so such a bit string oi . . . exists if and only if there is a c such that aj + e'^"^/^^ ^ (3j — (in fact 
c = {2^~^ + J2p=i ap2''~^)mod 2-'). Further, the inductive hypothesis ensures that for all / < j, 

= 0, 

thus such a bit string ai . . . exists if and only if \Cj\ 7^ 0. Hence, Cj is non-empty for j < ^ if and only if 
3ai . . . flfeVl < j < k{aj = e"^p=i "/^j). The condition that for j > k and aU ak+i . ■ .Uj G {0, 1}^"'' 



Uj 7^ e'^'5i;p=i '^p/^^ jg equivalent to \Cj\ — 0, since \Cj\ = requires that there exists a c such that 
aj + e^'"!"^''^ Pj = and -h e'^^'^/^'^Vfc 7^ 0. The only c < 2'= which satisfies this is c = Y%=i ap2P-\ so 
by symmetry any c which satisfies this must be able to be written as c = X]p=i ap2^'~^ for some Uk+i ■ ■ ■ aj. 
Hence we see that ^ and (pUjl are equivalent. □ 



Condition ([n]) of Lemma [T] corresponds to a more intuitive requirement for the transformed state to be 
well-formed. It says that for each j > 1 there must be a value of c such that /(c) = with the jth term in 
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([7]) equal to zero and the first j — 1 terms non-zero. If for some k there is no c satisfying this condition, then 
there must not be any c satisfying them for j > k either, or the state will not be well-formed. Condition 
(pH)) translates this notion into formal requirements about the relationships between the components of the 
untransformed input state components ai,f3i which will ensure the transformed state will satisfy condition 
(in]) and thus be well-formed. 

Lemma [1] gives us conditions for when the first condition of Theorem [T] is satisfied and it remains to 
determine which separable input states also satisfy the condition that VzHi^N)), with Z = X^fLi 2"^', is 
pair product invariant. The amplitudes which are deleted by the function Vz are the Z values of c which 
are in Cj for some j. 

Lemma 2. Let \iPn) be a separable input state for which the transformed state {i/jn) is well-formed, i.e. 



1=1 ■ 



The 



\iPn) satisfies the conditions of Lemma{^ Let k be as in Lemma{l\part ( fm)) 
2?z(|'0A')) is pair product invariant if and only if for all j > k + 1, aj/3j = 0, i.e. the (fc -f- l)th qubit can be 
in an arbitrary .superposition, and qubits k + 2 to n must not be in a superposition, although arbitrary phase 
is permitted. 

Proof. Let c' be the smallest c such that /(c) ^ 0, and let n' = n — k, N' ^2"\K = 2''. By symmetry, 
the N' non-zero amplitudes are f{dK + c') ioi < d < N' — 1. The zero-deleted state is thus VzHi^N)) — 
(/'(O), . . . , f'{N' — 1)), where f'{d) = f{dK + c'). By breaking up the sum in ([7]) we see that each of these 
amplitudes is of the form; 



1=1 



1=1 



(8) 



1=1 



where L 



2\ S 



c'/K is independent of d, as also is F = -^Yti^Mi + e^^""^^""" ' ^ Pi) ^ 
(recall k > I so dK/L is a positive integer). For all 1 < j < n', < mi < TO2 < >//2, pair product 
invariance (recall Definition [S]) requires that both f'{rni)f'{J — mi — 1) = f'{m2)f'{J — TO2 — 1) and 
/'(toi)/'(J/2 - mi - 1) = /'(m2)/'(J/2 - ma - 1). Since each f{d) ^ 0, we require 



/'( J - m2 - !)/'( J/2 - mi - 1) = /'( J - mi - l)/'(J/2 - ma - 1). 



(9) 



Symmetry means the left- and right-hand sides both contain common factors of F^, as well as j — 1 factors 
from the product for each transformed amplitude, due to the fact that ^^'^^J/^ — ^t^iJ/l -j-^j. i ^ j Thus 
the condition ^ simplifies to 

n' 
1=3 



1=3 



(10) 



which holds for all j, mi,m2 if and only if 'Dz{\''Pn)) is pair product invariant. 

We now show by induction that (ITOl) is satisfied if and only if for all 1 < j < n' , ak+j(3k+j — 0. Firstly, 
consider the case that j — n! . The products in (jlOp each contain only one factor, and expanding leaves only 
the cross-terms, and the condition simplifies to 



a. 



^l^^^^27ri{N'-m2)/N' _^ ^2^i{N' /2-7ni) /N' ^ Q;„^„ (e^''*'^^'""^ -|- q'^M^' /2-m2)/N' y ^j^-j^-j 
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Since this must hold for ah distinct mi,r7i2 only the trivial solution is possible, hence q;„/3„ = 0. 

Now, assume that ak+iPk+i = ior I — n' , . . . , j + 1, j > 1, and consider ak+j^ Pk+j- The products in 
pop run from j to n' , but all factors for I > j cancel when the pairs on each side are expanded since, by the 
inductive hypothesis, at+iPk+i = for these terms. The condition then reduces to a single factor and we 
find ak+jPk+j = exactly as in ([TT|). 

Hence, the transformed state is pair product invariant if and only if for all 1 < j < n' we have ak+j(3k+j = 
0. □ 

Theorem 2. Given a separable input state \4'n), the transformed state H^n) is separable if and only if 
(30 < fc < n)(3ai . . . e {0, 1}''') (vi <j<k aj ^ e"^'=i "'Z^'^' /3j 

A (afe+i ^ ±e"5:f=i -'/^'-'^'(ik+i) A (Vn > j > fc + 1) [a,l3, = 0]) . 
Proof. The proof follows directly from Lemmata [T] and [H □ 

Theorem [5] allows us to determine if a given separable state \iPn) will be entangled or not by the QFT. 
While the set of such states which are not entangled by the QFT is infinite, the conditions are still highly 
restrictive, and there is only one qubit that can ever truly be in an arbitrary superposition. However, the 
conditions between each ai and Pi are relative, so separability of the transformed state is invariant under 
phase rotations of individual qubits. These conditions, while restrictive, could be of value in developing new 
algorithms which make use of the QFT and give a strong insight into the entangling power of the QFT. 

4-3. Product-state De- quantisation 

For the set of states which are not entangled by the QFT, we can use the conditions of Theorem [5] 
to extend the basis-state de-quantisation. Let fc be as in Theorem m Let r = EjU a,.+, =o 2"^''+^'^ and 
uj — g27nr ^j^g coefficient of {ak+2 + Pk+2) ■ ■ ■ (ck,i + (in) in /(I). The de-quantised algorithm for states 
which are not entangled by the QFT is the following {b[x\ is the xth component of b starting from 0): 

Separable De-quantised QFT 

Input: The n two-component complex vectors ■ ■ - bn. 
Output: The n transformed vectors 6162 • • • bn- 

1. Calculate fc, ai . . .au as in Theorem [5] 

2. Calculate r, uj 

3. For j = 1 to fc -M: 

4. Set bn-jj^x = 

5. End For 

6. For j^lton — fc — 1: 

7. Let ; = n - j + 1 

9. End For 

10. For j = 1 to n: 

11. Set 6„_j+i[l] = wS„_j+i[l] 

12. Set w = 

13. End For 

Theorem 3. The Separable De-quantised QFT algorithm correctly computes the transformed n-qubit state 
\4'n) ~ Fn IV'Jv); where {iPn) is separable and the cth component of \'4'n) is described by ([7]), and does so in 
0{n) time. 
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Proof. The values of k and ai . . .Uk can be found readily in 0{n) time by just checking each pair aj, Pj 
too see which option, aj = 0,1, makes the first condition of Theorem [5] true, and setting aj accordingly, 
until neither is true, at which point k is found. Also, r and uj are efficiently found by direct calculation. It 
remains to verify that the algorithm correctly produces the state 



ft- 



1 

'7^ 



fe+1 



The algorithm calculates the amplitudes for each qubit, so if we let the n-bit binary expansion of c be 
c„ . . . ci we have 



/(c) = 6i[c„] • 62[c„-l] • • • bn[ci] 



'k+i 



Note that, since a,- = or /3, = 0, 



n K-+ft) 

j=k+2 



n K+ft) 

]=k+2 



(12) 



n( 

j=k+2 



clV-\ 



n (a,+/3,)=w^ n K-+ft), 

j=fe+2 J=fe+2 



so our algorithm produces this term correctly. 

Since the output state is separable, the conditions of Theorem [5] must be satisfied and only one out of 
the first K amplitudes is non-zero. This amplitude is the one with c' = X^zLi ^(2'"^, and by symmetry all 
the other non-zero amplitudes occur at c = c' + d2^^^ for < d < iiT — 1. To verify this, note that for all 
j < k, we have 



I2i'=i ai2 



ft ^ 0, 



and hence /(c') ^ 0. From it is clear that /(c) is calculated correctly for these values of c. For all other 
values of c which have ci . . . c„ 7^ ai . . . a„, let m be the smallest i < n such that Ci ^ ai. Then we have 



ftn — 



/3m =0, 



and hence /(c) is correctly produced for all c. 

The algorithm is also clearly seen to require 0{n) time, and thus the proof is completed. 



□ 



This algorithm has all the advantages of the basis-state de-quantised algorithm, but operates on a much 
larger ranger of input states, making it a much more powerful de-quantisation. Importantly, just like the 
basis-state de-quantisation, it is actually more efficient than the QFT algorithm. While this algorithm 
will not work on all separable input states like the tensor-contraction simulation in Q, it is a stronger 
de-quantisation in the sense that it gives a complete description of the output state as opposed to the 
probability of measuring a particular value, and is trivial to use as a subroutine in a larger de-quantisation. 



5. Discussion 

The ability to de-quantise the QFT algorithm brings up some interesting points. The two de-quantisations 
presented in this paper compute the Fourier transform on a restricted set of input states. On the other hand 
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the standard QFT algorithm computes the Fourier transform on arbitrary separable or entangled input 
states. In fact, the standard QFT algorithm is a quantum implementation of the basis-state algorithm, 
but the linearity of quantum mechanics ensures that arbitrary input states are transformed by this simple 
algorithm. De-quantisation techniques such as the one presented, as well as those of 0, 0, 0|i have to 
efficiently simulate the linearity that is inherit in the quantum mechanical medium. The de-quantisations in 
this paper highlight the important distinction that should be made between the quantum Fourier transform 
and the quantum algorithm computing it. The QFT is a unitary transformation of an n-qubit state, while 
the QFT algorithm is a recipe for creating a sequence of local gates which computes the QFT on a given 
state. While these two notions are equivalent in quantum computation, when we depart from quantum me- 
chanics this is no longer the case, and the de-quantised algorithm does not suffice to compute the complete 
QFT. 

It is interesting to note that both de-quantisations presented in this paper run in 0(n) time, more 
efficient than the 0{n^) of the quantum algorithm. This is due to the restrictions imposed by measurement 
no longer being present when we develop a classical counterpart. This increase in efficiency is something 
not seen in other de-quantisations of the QFT which are based on the quantum circuit topology, and thus 
inherently and perhaps unnecessarily work within the restrictions the quantum circuit was designed under. 
The Separable De-quantised QFT algorithm computes the QFT on a large number of input states and it 
remains to be seen if this de-quantised algorithm can be applied to existing or new quantum algorithms to 
produce further de-quantisations. The fact that both the input and output states are separable also ensures 
the existence of a de-quantised inverse algorithm too, which is of practical significance. 

Another issue worth noting is that we must be careful to consider the complexity involved in manipulating 
the complex amplitudes in a state-vector when performing de-quantisation. While it did not contribute to 
the complexity of the de-quantised algorithms presented in this paper, attention had to be paid to make 
sure this was the case, as this would not have been so if we implemented the directly obvious algorithm. In 
quantum computation, however, the amplitudes are just our representation of a property of physical states. 
It is these physical states, rather than the amplitudes, which are altered by unitary transformations, and as 
a result we observe the amplitudes changing. This reiterates the need for care when de-quantising, as the 
amplitudes have no a priori reason to be easily calculated, or computable at all for that matter. 

6. Summary 

We have shown that the quantum algorithm computing the QFT can be de-quantised into a classical 
algorithm which is more efficient and in many senses simpler than the quantum algorithm, primarily because 
the need to avoid measurement of the system is no longer present. However, the direct de-quantisation of 
the QFT algorithm leads to a classical algorithm which only acts on a basis-state. This difference is due 
to the linearity of quantum computation ensuring a basis-state algorithm computes the complete QFT, 
highlighting this linearity as a key feature in the power of the QFT. By examining the entangling power of 
the QFT we devised conditions for when the QFT leaves a separable state unentangled, and showed that 
this separability is invariant under phase-rotation of the input qubits. We extended our de-quantisation to 
work on this set of states without any loss of efficiency. 

This de-quantisation of the QFT serves not only to illustrate more deeply the nature of the QFT, but also 
provides a useful tool for possibly de-quantising algorithms using the QFT with very little effort. Further, 
the techniques involved can help identify de-quantisable algorithms more easily, as well as aiding the creation 
of new quantum algorithms and subroutines by deepening our understanding of what is needed to make a 
quantum algorithm so useful. 
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